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A class of exact solutions is obtained for the Lienard type ordinary non-linear differential equation. 
As a first step in our study the second order Lienard type equation is transformed into a second 
kind Abel type first order differential equation. With the use of an exact integrability condition 
^ ■ for the Abel equation (the Chiellini lemma), the exact general solution of the Abel equation can be 

obtained, thus leading to a class of exact solutions of the Lienard equation, expressed in a parametric 
form. We also extend the Chiellini integrability condition to the case of the general Abel equation. 
As an application of the integrability condition the exact solutions of some particular Lienard type 
equations, including a generalized van der Pol type equation, are explicitly obtained. 
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I. INTRODUCTION 

The Lienard type second order nonlinear differential equation of the form [l| 

<n : 

x{t)+f(x)x{t)+g{x) = Q, (1) 

\Q ; 

where f(x) and g(x) are arbitrary real functions of x, with f(x),g(x) £ C°°(J), defined on a real interval I C 5fi, as 
well as its generalization, the Levinson-Smith type equation Q 

<N i x{t) + f(x,x)x(t)+g(x) = 0, (2) 

o ■ 

■ where a dot represents the derivative with respect to the time t, and / is a function of x and x, plays an important 
role in many areas of physics, biology and engineering The mathematical properties of these types of equations 
have been intensively investigated from both mathematical and physical points of view, and their study remains an 
• active field of research in mathematical physics (34l3j . 

From a physical point of view, the Lienard equation represents the generalization of the equation of damped 



oscillations, 'x + jx + lu 2 x — 0, where 7 and uj 2 are constant parameters, respectively 14|. For 7 = we obtain the 
equation of the linear harmonic oscillator, which represents one of the fundamental equations of both classical and 
quantum physics. Generally, a linear oscillation can be described by the equation x + f(t)x + g(t)x = 0. A particular 
type of the general Lienard equation, the van der Pol equation [l5[ 

x(t)- fj,[l- x 2 (t)]x(t)+x(t) =0, (3) 

describing a non-conservative oscillator with non-linear damping, is extensively applied in both the physical and 
biological sciences. The Lienard type equations can also be used to model fluid mechanical phenomena. The linearly 
forced isotropic turbulence can be described in terms of a cubic Lienard equation with linear damping of the form [16( 

x(t) + [ax(t) + b] x(t) + cx(t) - x 3 (t) + d = 0, a, b, c, d = constant. (4) 
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Lienard type equations also naturally appear in the mathematical description of some important astrophysical 
phenomena. For example, the time-dependence <j>(t) of the perturbations of the stationary solutions of spherically 
symmetric accretion processes can be described by a generalized Lienard type equation of the form [13] 

t + eff^J^ + v'^) = o, (5) 

where e is a small parameter, and V(<p) is the potential of the system, with the prime indicating the derivative with 
respect to (f>. A dynamical systems analysis of this Lienard equation reveals a saddle point in real time, with the 
implication that when the perturbation is extended into the nonlinear regime, instabilities will develop in the accreting 
system. 

It is the purpose of the present paper to introduce some exactly integrable classes of the Lienard equation, Eq. ([IJ , 
whose solutions can be obtained in an exact analytical form, and to formulate the integrability condition for this class 
of differential equations. To obtain the functional form of the integrable Lienard type equations we reduce them first 
to an Abel type equation (l8| . Then, we apply to the latter Abel equation an integrability condition, equivalent to 
the initial Lienard equation, that was obtained by Chiellini [ljj[2(|. In fact, the Chiellini condition has been recently 
used for the study of the Abel differential equations, as well as for the second order differential equations reducible 



to an Abel type equation, in [2l[- 25[. In [24| the differential Chiellini integrability condition was reformulated in an 
integral form, and the general form of the solution of the Abel equation was obtained in a simpler form. 

Therefore it turns out that, if the coefficients f(x) and g(x) of the Lienard equation satisfy two specific conditions, 
then the general solution of the Lienard equation can be obtained in an exact analytical form. Some examples of 
exactly integrable Lienard equations, of physical interest, are also considered. The generalization of the method to 
the case of the Levinson-Smith type equations of the form ^ is briefly discussed. 

The present paper is organized as follows. In Section [Til we introduce the Abel equation representation for the 
Lienard equation, and we formulate the integrability condition of the first order Abel equation. In Section Mil we 
obtain the general solution of the Lienard type equations satisfying the integrability condition of the Abel equation. 
The exact solutions of some non- linear Lienard type differential equations are obtained in Section HVl We discuss and 
conclude our results in Section fVl 



II. REDUCTION OF THE LIENARD EQUATION TO AN INTEGRABLE ABEL TYPE EQUATION 

As a first step in our study of the Lienard equation (fTJ) we reduce it to an Abel type first order non-linear differential 
equation. Then, by using an integrability condition for this equation, which involves a differential relation between 
the coefficients f(x) and g(x) of the equations, we obtain the general solution of the Lienard equation in an exact 
parametric form. 

By denoting x = u, the Lienard equation (flj can be written as 

du 

u— + f(x)u + g(x)=0. (6) 

By introducing a new dependent variable v — 1/u, Eq. ([6]) takes the form of the standard second kind Abel differential 
equation, 

^L = f( x y+g( x y. (?) 

ax 

In this context, an exact integrability condition for the Abel equation Eq. ([7]) was obtained by Chiellini [l9| (see 
also and can be formulated as the Chiellini Lemma as follows: 

Chiellini Lemma. If the coefficients f(x) and g(x) of a second kind Abel type differential equation of the form 

^ = f( x y + g ( x y 7 ( 8 ) 

ax 



satisfy the condition 

d g(x) 



dx /(a 



kf(x), k — constant ^ 0, (9) 



then the Abel Eq. (JSJ can be exactly integrated. 
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In order to prove the Chiellini Lemma we introduce a new dependent variable w defined as 19, 2 



v = —r^rW. 
9\x) 



Then Eq. JS} can be written as 

1 df(x) fjx) dg(x) 
g(x) dx 9 2 {x) dx 

On the other hand the condition given by Eq. Q can be written in an equivalent form as 

f{x) dg(x) 1 df(x) _ k f(x) 



g(x) ax g (x) 



Therefore Eq. (jlip becomes 



g 2 (x) dx g[x) dx 
dw P(x) 



g 2 {x) 



(w 2 + w + k) , 



dx g(x) 

which is a first order separable differential equation, with the general solution given by 



dx = 



dw 



w (w 2 + w + k) 



= F(w,k). 



With the use of the condition © the left hand side of Eq. (ITO can be written as 24 1 



f 2 ( x ) dx= 1 
g(x) k 



dx 



hi 



g(x) 



dx = — In 
k 



/(a 



Co, 



where Co is an arbitrary constant of integration. Therefore the general solution of Eq. (1131) is obtained as 

g(x) 



/(*) 



= C _1 e F(tu ' fc) , 



where C 1 = exp (— fcCo) is an arbitrary constant of integration, and 



\/w' 2 -\-w-\-k 



exp 



i 



l+2w 



2arctanh(l + 4w) 



k = 



+2m 



-l/2Vl-4fe 



l+2tu 
Vl-4fe 



1/2V1-4S; 



fc < 1/4, 



respectively. Eq. (| determines w as a function of x. 

The integrability condition given by Eq. (|9]) can be written as 



dgjx) _ 1 df(x) 
dx f(x) dx 



g(x) + kf 2 (x), 



(10) 



(11) 



(12) 



(13) 



(14) 



(15) 



(16) 



(17) 



(18) 



representing a first order linear differential equation in g. As a function of g, the function / satisfies the differential 
equation 



1 df(x) 



1 dg(x) 



f(x) dx g(x) g(x) dx 

In order to solve Eq. (fH?)) we introduce a new dependent variable fix) = l/a(x). Therefore Eq. (|19p becomes 

k 



da(x) 1 dg(x) 2 
a \ x ) 1Z~ = 173 j. g W 



g(x) dx 



(19) 



(20) 
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By denoting <r 2 (x) = £(x), we obtain a first order differential equation for £, 



dx 



g(x) dx 



2k 

9{x)' 



(21) 



Therefore the Chiellini Lemma can be reformulated as: 

Lemma 1. If the coefficients f(x) and g(x) of the Abel Eq. (J5J satisfy the conditions 



or 



g(x) = f(x) 



m = ±- 



d+k f(x)dx 



(22) 



(23) 



■\Jc 2 + 2kJ g(x)dx 

where C\, C 2 , and k are arbitrary constants, the Abel equation is exactly integrable, and its solution is given by 

v(x) = Ce- Fiw ^^w(x), (24) 



where the functions F(w, k) are given by Eqs. (fl7l). A similar result was obtained in [24 
The Chiellini Lemma can be extended to the general Abel equation of the form 



dv 
dx 



a(x) + b(x)v + f(x)v 2 + g(x)v 3 



(25) 



where a(x), b(x), f(x) 7 g{x) G C°°(I) arc defined on a real interval ICS, and a(x),b(x) ^ 0,Vx € /, as follows. By 
introducing a new function p(x), defined as 



Eq. becomes 



v(x) = eS b{x)dx p{x), 



= a(x)e- S b ^ dx + f(x)eJ b ^ dx p 2 + g(x)e 2 / 6 ^ V- 

ax 



We assume now that the functions 6(x), /(x) and <?(x) satisfy the condition 

d 9 {x)ef»W* _ fb(x)dx 
dx f(x) ~ Kl ^ x > e 

where k\ is an arbitrary constant. Then, by introducing the transformation 



p(x) 



g(x) e Sb{x)dx 



s(x), 



Eq. ((27]) becomes 



ds g(x) , . f 2 (x) , , \ 
' v \a(x) + L j-^-s (s 2 + s + fci) . 



dx /(x) 



(26) 



(27) 



(28) 



(29) 



(30) 



Hence we have obtained the following generalization of the Chiellini Lemma: 

Lemma 2. If the coefficients of the general Abel Eq. (f25j) satisfy the conditions (|28|) and 

/ 3 M 



a(x) = h 



g 2 {x) 



(31) 



where fc 2 is an arbitrary constant, then the Abel equation can be exactly integrated, and its general solution is given 
I'.v 



( \ ^( x ) ( \ 
V{X) = -gTxj^ 



(32) 
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with s(x) a solution of the equation 

G (s, k-L,k 2 ) 



dx = — In 



g(x) h 

where Kq is an arbitrary constant of integration, and 

G (s, fci,fc 2 ) = 

By using Eq. (gHJ), Eq. fl33]) becomes 

g(x) e J b ^ d ' 



g(x)ef b ^ dx 



/(*) 



ds 



s 3 + s 2 + k\s + A?2 



/(*) 



K , 



(33) 



(34) 



(35) 



where A"i = exp (— kiK$) is an arbitrary constant of integration, and G (s, fci, ^2) = ^iGq (s, fci, fe), respectively. 



III. A CLASS OF EXACT SOLUTIONS OF THE LIENARD EQUATION 

As we have already seen, the second order non-linear Lienard Eq. (TTJ) can be reduced to an Abel type equation of the 
form given by Eq. with the general solution given by v(x) — C exp [— F(w(x), k)] w(x), where w(x) is determined, 
as a function of x, by Eq. (|16p . Alternatively, Eq. (fT6|) fixes x as a function of w, 

(36) 



x = x(w) 

To find the time dependence of x, we start from 

dx dx dw 



1 g(x) 1 



which gives 



dt dw dt v f(x)w' 

dw dw g(x) 1 



dt dx f{x) w 

With the use of Eq. (p~3|) . satisfied by the function w(x), we obtain for dw/dt the equivalent expression, 

= f{x) (w 2 +w + k). 

Therefore we have obtained the following: 

Theorem. If the coefficients of the Lienard equation ([1} satisfy the conditions 



g(x) = f(x) 



or 



(37) 



(38) 



(39) 



m - ± 



Gi +k J f{x)dx 
^Jc 2 + 2kJ g(x)dx 



(40) 



(41) 



where C\, C2 and k are arbitrary constants, then the general solution of the Lienard equation Eq. (fl} can be obtained 
in an exact parametric form, with w taken as parameter, as 



dw 



f(x(w)) (w 2 + w + fc) ' 
with x(w) obtained as a solution of the equation 



x = x(w), 



9{x) 
fix) 



c -l e F{ w ,k) ; 



(42) 



(43) 
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and 

F(w,k) = k [ ^ (44) 

J w (w z + W + K) 

Similar results were previously obtained in [24| , where the integral form of the Chiellini integrability condition were 
explicitly formulated. A particular integrable case of the Lienard equation can be obtained for the case k = 0. In this 
case, the Chiellini condition given by Eq. ^ immediately provides 

g(x)=Af(x), (45) 

with A an arbitrary constant. Thus, the Lienard equation takes the particular form 

x + f(x)x + Af(x) = 0, (46) 

with the associated Abel equation given by 

^ = f( X y(l + Av), (47) 
where x — 1/v. The general solution of Eq. (j47|) is given by 

f(x)dx = Aln(-+A)--+Ki, (48) 
\v J v 

where K\ is an arbitrary constant of integration. Therefore, the general solution of Eq. (1461) can be written in a 
parametric form, with v taken as a parameter, in the following form 

t-to = f f, i vTTi , AV x = x{v), (49) 
J f{x(v))v{l + Av) 



where x = x(v) is the solution of Eq. 

In the general solution for the time, given by Eqs. (|42|) and (|49)) . one can take the arbitrary integration constant 
to as zero, without any loss of generality. This choice fixes the origin of time at t = 0. The arbitrary integration 
constant C, as well as the initial value wq of the parameter w can be determined from the initial conditions at t = 
for the position x and the velocity x, given by 

a:(0)=a:o, i(0) = x , (50) 
where xq and Xq are the initial values of x and x at t — 0. By evaluating Eq. (|43[) for x — Xq, we obtain 



while evaluating Eq. (|3T[) at t = gives the equation xq — [g (xq) / f (xo)] w 1 , which determines the initial value of 
the parameter wq as 

«*=-MH (52) 

Once the initial value of the parameter wq is known, the value of the integration constant C^ 1 is obtained as 

c- 1 = i^A e -F( WOtk )_ (53) 

The procedure for the exact integration of the Lienard type equations based on the Chiellini Lemma can be easily 
extended to the generalized Lienard equations of the Levinson-Smith form, given by Eq. ([2]), if they can be transformed 
to an Abel type equation. As a particular case of the integrable Levinson-Smith type non-linear differential equations 
we consider the equation 



x + [-f(x)x 2 + S(x)x + f(x) \ x + g(x) = 0, 



(54) 
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where j(x) and S(x) are some arbitrary functions of the variable x. By denoting x — 1/v, Eq. (|54|) takes the form of 
the general Abel equation 

dv 

— = ■y(x) + S(x)v + f(x)v 2 + g(x)v 3 = 0. (55) 

If ^(x) — 0, by means of the transformation v{x) = e$ s ^ dx h(x), Eq. (|55j) can be written in the standard form of 
the Abel equation, 

^ = A{x)h 2 + B{x)h 3 , (56) 

where A{x) = f(x)e^ s ( x ) dx } and B(x) — g{x)e 2 $ s ( x ) dx . If the coefficients A(x) and B{x) of the equation satisfy the 
conditions of Lemma 1, then the general solution of Eq. (|56l) can be obtained through quadratures. If ^(x) ^ 0, then 
from Lemma 2 it follows that if the functions ^(x), 5(x), f(x) and g(x) do satisfy the conditions 

d n(r)pf s ( x ) dx f 3 (x) 

-r 9 [ ' ' > = k 1 f(x)ef 5 ^ dx , 7 (*) = fc 2 LW h,k 2 = constant, (57) 

dx f(x) g 2 {x) 

then the generalized Lienard type equation (|54p can be integrated exactly. Therefore all the intcgrability results 
obtained for the Lienard equation can be applied for the Levinson-Smith type equations of the form (|54[) . 

IV. EXAMPLES OF EXACTLY INTEGRABLE LIENARD TYPE EQUATIONS 

In the present Section, we consider some exactly integrable Lienard type equations, which represent the general- 
izations of Eqs. © and dU, respectively. As a first case we assume that the functional form of the function f(x) is 
known. Then the Chiellini intcgrability condition fixes the form of the function g(x), and allows to find the general 
solution of the Lienard equation in an exact parametric form. The case in which the function g(x) is fixed is also 
considered. Furthermore, an integrable generalization of the van der Pol equation is also explored. 

A. First case: f(x) = ax + b. 

As a first case we assume that the function f(x) is given by 

f(x) = ax + b, a,b = constant. (58) 
Then from the first integrability condition, given by Eq. (|40[) we obtain the function g{x) as 

1 3 

g [x) = -a 2 kx 3 + -abkx 2 + (aCi + b 2 k) x + bCx, (59) 

where C\ and k are arbitrary integration constants. Therefore the exactly integrable Lienard equation is given by 

1 3 

x + (ax + b) x + -a 2 kx 3 + -abkx 2 + (ad + b 2 k) x + bCi = 0. (60) 
As a function of the parameter w, x is determined by Eq. (|43|) . which gives for x the quadratic algebraic equation 



^x 2 +bkx + C\=C- 1 e F ( w > k \ (61) 



which determines function of w as 



x(w) = v - ; . (62) 

ak 



-bk ± ^b 2 k 2 - 2ak [d - C- 1 e F ( t °. fe ) 
a 

The time dependence of x is determined as a function of w as 

/kdw ,„„, 
; = ■ (63) 

Jb 2 k 2 - 2ak [d - C- 1 e F ( t °. fe )] (w 2 + w + k) 

Eqs. (|6"2")l and (|53"]) give the exact solution of the Lienard Eq. (|60p. Depending on the value of the constant k there 
are three distinct classes of solutions of this equation. 



B. Second case: g(x) — cx + d 

Secondly, we consider the case in which the function g(x) is fixed. By analogy with the van der Pol Eq. ([3]), we 
assume that 

g{x) = cx + d, c,d — constant. (64) 
Then, after determining the function f{x) from the integrability condition Eq. (1411) . we obtain the Lienard equation 

cx ~\~ d 

x ± : x + cx + d — 0. (65) 

y/ckx 2 + 2dkx + C 2 



Eq. (j43|) gives the equation 



with the solution 



ckx 2 + 2dkx + C 2 =C- 2 e 2F{w > k \ (66) 



-dk ± Jd 2 k 2 - ck \C 2 - C- 2 e 2ir O> fc ) 

x(w) = 1 -. (67) 

ck 

The parametric time dependence of the solution is obtained as 

k f e F{w ^dw , s 

t-t a = ±- . (68) 

° J ^d 2 k 2 - ck [C 2 - C- 2 e 2i? K*0] {w 2 +w + k) 
Eqs. (|67| and (|68)l give the exact analytic solution, in a parametric form, of the Lienard Eq. (|65| . 



C. Third case: the generalization of the van der Pol equation 

Finally, we consider the integrable generalization of the van der Pol Eq. ([3]), in which we fix the function f(x) as 
f(x) = —fi (l — x 2 ), and obtain the function g{x) from the integrability condition Eq. (|40[) . Therefore the integrable 
generalization of the van der Pol equation is given by 

x - [i (l - x 2 ) x + \k^ 2 x 5 - ^kfi 2 x 3 + dfix 2 + kfi 2 x - Cin = 0. (69) 



The parametric dependence of x is determined from the algebraic equation 

Cx - kpLX + h^x 3 = C- 1 e F{w ' k) . 



(70) 



Eq. (|70p can be rewritten in the form 

x 3 — 3x + H (F) = 0. (71) 

where we have denoted 



The solution of the algebraic Eq. (JTTJ) are given by 



2 i/3 { V^ 2 [F(w, k)] - 4 - H [F(w, k)]j 



1/3 



x(w) = — — - + . (73) 

| y^H 2 [F(w, k)]-4-H [F(w, k)} \ 

In order to have a real solution of the cubic Eq. (|7ip . the conditions 

H 2 [F(w, k)} - 4 > 0, (74) 
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and 

^H 2 [F(w, k)]-A-H [F(w, k)} > 0, (75) 

must be satisfied for all w. 

The parametric time dependence of the solution of the generalized van der Pol equation is obtained as 

1 f 4> 2 {w,k)dw 
1 ~ *° = m J [^{w, k) + <f> 2 (w, k) + 1] (w 2 + w + k) ' (76) 

where we have denoted 

2 l/3 

[w, k) = — -rj^- (77) 

^/H 2 [F(w,k)] -4 - H [F(w, k)} 

Depending on the numerical values of the parameters k, /i, C\ , C a large class of dynamical evolutions of the solutions 
of the generalized van der Pol equation can be obtained. 

V. DISCUSSIONS AND FINAL REMARKS 

In the present paper we have introduced a class of exactly integrable Lienard, and generalized Lienard type 
equations. If the coefficients of the second order non-linear equations satisfy some specific conditions, which follow 
from the Chiellini Lemma, then the general solution of the Lienard differential equation can be obtained in an 
exact parametric form. As an application of the integrability procedure obtained, we have considered some specific 
examples of exactly integrable non-linear differential equations that could be of physical interest. One of these 
equations, Eq. (|60|) is similar in form with Eq. (|4]), and in fact represents the exactly solvable generalization of 
the equation describing the linearly forced isotropic turbulence [![. We have also considered an exactly solvable 
generalization of the classical van der Pol oscillator equation, in which higher order force terms are also present. In 
all these cases of physical interest the general solution of the corresponding Lienard equation can be obtained in 
an exact parametric form. The existence of an analytical solution may allow a deeper understanding of the highly 
non-linear physical processes that govern most of the natural phenomena. 

The exact solutions also allow us to obtain some approximate solutions of the considered differential equations, 
corresponding to the small and large values of the parameter w, respectively. In the limit of small w, i.e., w <C k, 
giving exp [F(w, k)] s» w, Eq. (j4"3)) takes the simple form 



g(g) 

while the parametric time evolution can be obtained as 



CrV (78) 



1 f dw 

'-^kllWTY (79) 

In the limit of large w, so that w ^> k, and w 2 w, exp [F(w, k)} ps exp (k J dw/w 3 ) = exp \—k/2w 2 ), and the 
approximate asymptotic solution of the exactly integrable Lienard equation is given by 



and 



C- l e- k/2w \ (80) 



/dw 
ff ( » 2 - (81) 

As an application of the previous asymptotic equations we consider the case fix) = ax + b, with a,b = constant, 
jiving g(x)/f(x) = C\ + kax 2 /2 + kbx. 
In the limit of small x, by neglecting the x 2 term, we obtain 
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and 



6C, 
t — Iq ~ — In 
a 



aw 



ad + o 2 fc 



respectively, giving 



x(t) 



o a(t-t )/bC 



abk a 

In the limit of large x, so that kbx 3> C%, and a:r/2 ^> &, respectively, we obtain g(x)/ f(x) 



kax 2 /2, and 



x(w) 



2C- 1 
ka 



-k/i-b 



(83) 



(84) 



(85) 



t-t 



dw 



2a/kC e ~ k / iw2 + b 



(86) 



In the range of the values of w for which ^/2a/kC i 



-/c/4u> 2 



3> 6, we obtain 



lirC 



Vk^ 



t - t m - \ — erfi — , 



2a 



2w 



(87) 



where erfi(,z) gives the imaginary error function erfi(z) = erf(iz)/i. In the large time limit the solution of the Lienard 
Eq. (|60[) can be obtained only in a parametric form. 
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